IFT-UAM/CSIC-04-44 
CI-UAN/04-07FT 
hep-ph/0409070 



2 



Symmetry Nonrestoration at High 
Temperature in Little Higgs Models 

J.R. Espinosa^* M. Losada^^and A. Riotto^* 

^ IFT-UAM/CSIC, Fac. Ciencias UAM, Cantoblanco, 28049 Madrid, Spain 
Centro de Investigaciones, Univ. Antonio Narino, Cll. 58A No 37-94, Bogota, Colombia 

Q 

INFN, Sezione di Padova, via Marzolo 8, Padova 1-35131, Italy 

February 2, 2008 



Abstract 

A detailed study of the high temperature dynamics of the scalar sector of Lit- 
tle Higgs scenarios, proposed to stabilize the electroweak scale, shows that the 
electroweak gauge symmetry remains broken even at temperatures much larger 
than the electroweak scale. Although we give explicit results for a particular 
modification of the Littlest Higgs model, we expect that the main features are 
generic. As a spin-off, we introduce a novel way of dealing with scalar fluctua- 
tions in nonlinear sigma models, which might be of interest for phenomenological 
applications. 
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1 Introduction 



Little Higgs (LH) models [1-4] provide a new scenario for electroweak symmetry break- 
ing which stabilizes the electroweak scale. It is well known that symmetries in the 
Standard Model (SM) Lagrangian force gauge bosons and fermions to be massless. 
The main feature of LH models is the inclusion of additional global symmetries which 
force the Higgs boson mass to be zero. The incorporation of explicit violations of these 
symmetries in a precise way explains why the ratio of the Higgs mass to the cutoff of the 
theory, m|^/A^, is small: the Little Higgs boson is a naturally light pseudo-Goldstone 
boson [2]. 

Below the cutoff scale A ~ 10 TeV, there are additional particles (with masses of 
order / ~ 1 TeV) which cancel the quadratically divergent contributions to the Higgs 
boson mass from SM particles (or at least the most dangerous ones) . Beneath the TeV 
scale the effective degrees of freedom are those of the Standard Model. The initial 
implementation of this idea was based on a SU(5)/SO(5) nonlinear sigma model which 
contained a gauged [SU{2) x f/(l)]^ subgroup. However, this model is in conflict with 
low-energy precision electroweak measurements and from direct searches for a Z' boson 
[5-9], both problems being related to the additional U{1) gauge group. Although this 
problem is not of direct concern for our purposes, we focus in this paper on a Little 
Higgs model which is a simple variation of the so-called Littlest one, having only a 
gauged SU{2) x SU{2) x U{1) group [10]. The performance of the model with respect 
to electroweak fits is improved while the smallness of the g' coupling tames the remnant 
quadratic divergence in the Higgs boson mass associated to the U{1)y gauge group. ^ 

In this paper we study the finite temperature behaviour of Little Higgs models. 
The main motivation is related to the peculiar way in which such models cancel the 
quadratically divergent corrections to the Higgs boson mass. As is well known, in the 
early Universe the Higgs scalar flcld gets an effective thermal mass (squared) trilj^jr ~ 

from interactions with the ambient hot plasma. It is this positive m^^j which is 

responsible for symmetry restoration at high temperature [12]. This effective mass can 

be computed diagrammatically using well established finite temperature techniques; 

-'^An interesting alternative has been proposed recently [11] in which a T-parity is imposed on the 
nonlinear sigma model such that it eliminates the strongest constraints from tree-level processes to 
electroweak observables. In any case, most of our findings are expected to be generic in Little Higgs 
models and should hold also for models with T-parity. 
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the self-energy diagrams which are quadratically divergent at T = are precisely the 
ones that give a contribution to the Higgs thermal mass. The correspondence is very 
simple [13]: a one-loop bosonic self-energy diagram that gives 5m1 = /tA^/(167r^), 
where A is an UV cutoff, produces 5^e// — 't^^/lS at finite temperature. For a 
fermionic one-loop self-energy diagram, if the zero T result is Sm\ ~ —kIs?/ (IGtt^), the 
finite T contribution of the same diagram will be Smlff — For instance, in 

the Standard Model one gets the well known quadratically divergent correction to the 
Higgs mass 

5?A^ r 1 , 1 _ „ _i 

(1) 



5ml 



167r2 

where A is the Higgs quartic coupling (normalized so as to have ml = 2Af ^ with v = 246 
GeV), g is the SU{2) gauge coupling, g' the U{1)y coupling and ht the top Yukawa 
couphng. The result (1) translates, at T 0, into 



(2) 



In the MSSM, although the sum of quadratically divergent corrections to ml cancels 
between bosonic and fermionic degrees of freedom, the corresponding corrections at 
T 7^ do not and one gets symmetry restoration also in that case. 

The defining property of Little Higgs models is that quadratically divergent cor- 
rections to ml cancel sector by sector in the model, between particles of the same 
statistics. At T ^ this leads one to expect that 5mlff will also be zero in this type of 
models. Notice though that this will only happen for sufficiently large temperatures, 
such that the heavy partners of the SM particles are thermally produced and populate 
the hot plasma. At low temperatures we expect the heavy particles with mass ~ / 
to be Boltzmann decoupled, and we find the same SM dependence on the Higgs 
mass, with the electroweak symmetry being restored as usual. (That is, the electroweak 
phase transition will take place as in the SM, at T ~ 100 GeV.) As the temperature 
approaches / ~ 1 TeV the new particles introduced in Little Higgs models will be in 
thermal equilibrium in the plasma, the thermal Higgs mass will drop to zero and the 
electroweak symmetry will be broken again^. Although these expectations will prove 
to be correct, they concern the behaviour of the Higgs potential, V{h), for small values 

^It is interesting to confront this result with the negative expectations of ref. [14], which claimed 
that symmetry nonrestoration could only be obtained at the price of hierarchy problems, while we 
find symmetry nonrestoration precisely due to the good UV properties of Little Higgs models. This 
contradiction is resolved by noting that the results of [14] do not apply in nonlinear sigma models. 
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of the Higgs field h. In order to study what happens to the minimum of the Higgs 
potential at T ~ / ~ 1 TeV, we will need to study larger values of h. 

Without knowing the UV theory that supersedes the Little Higgs model beyond 
A ~ 47r/ we cannot study the behaviour of our theory at T > A/tt ~ 4/ (at such 
temperatures the free-energy contribution from particles with mass ~ A is no longer 
negligible). With this limitation in mind, we would like to explore in this paper the 
behaviour of the Higgs potential of Little Higgs models at finite temperature. We will 
confirm the expected behaviour described in the previous paragraph and discover some 
pecuhar features in the temperature evolution of V{h). We present detailed results for 
the model of ref. [10] but expect that the main features we find are generic as they are 
based on the defining properties of Little Higgs models. 

In section 2 we present the model and our notation. Section 3 describes in de- 
tail the structure of the zero temperature effective potential at one-loop order, paying 
particular attention to the treatment of the contributions from scalar fluctuations. Sec- 
tion 4 then goes on to compute the finite temperature effective potential and describes 
some interesting features of its temperature evolution. In particular, we discover that 
the electroweak gauge symmetry remains broken even if the system is heated up to 
temperatures larger than the electroweak scale. In section 5 we present some conclu- 
sions. Appendix A contains explicit expressions for the mass matrices of the different 
species of particles in the model, calculated to all orders in the Higgs background, 
which are necessary to compute the one-loop effective potential. Appendix B deals 
with the behaviour of the potential along the direction of the triplet field contained in 
the model. 

2 The Little Higgs Model 

The model is based on the SU{5)/ S0{5) nonlinear sigma model of [2] (the Littlest 
Higgs), modified according to [10]. The spontaneous breaking of SU (5) down to 5*0(5) 
is produced by the vacuum expectation value (vev) of a 5 x 5 symmetric matrix $ 
[which transforms under SU{5) as $ — > C/$C/^] for instance when ($) = I5 (we denote 
by /„ the nxn identity matrix). This breaking of the global SU (5) symmetry produces 
14 Goldstone bosons among which lives the scalar Higgs field. Instead of working on 
the background ($) = I5 we follow [2] and, making a basis change, we choose ($) = Eq 



3 



where 

/ /2 \ 

Eo= 1 . (3) 
V ^2 J 

Calling Uo the SU (5) matrix that performs this change of basis, we have Sq = UqUq 
while all the group generators change as ta — Uot^^^Uo [4°-* are the generators in the 
original basis]. The unbroken SO (5) generators satisfied the obvious relation T^°^ + 
T^^^'^ — and multiplying on the left by Uq and on the right by Uq we arrive at the 
condition 

r„Eo + SoTj = , (4) 

for the generators in the new basis (a condition which is immediate to obtain alterna- 
tively just by requiring invariance of Eq). In the original basis the broken generators 
obviously satisfy = Multiplying again by Uq and Uq one gets in the 

transformed basis 

Xa^O — ^oXj ■ (5) 

The Goldstone bosons can be parametrized through the nonlinear sigma model field 



(6) 



where U^J^a n"X". The model assumes a gauged SU{2)i x SU{2)2 x U{1)y subgroup 
of S'C/(5) with generators 

a"" /2 02x3 \ r^a ( OsxS O3X2 



^'~[ 03X2 03x3 j ' I 02x3 '^'"'/^ 



(7) 



(where cr" are the Pauli matrices) and Y — diag(l, 1,0,— 1,— 1)/2. The vacuum expec- 
tation value in eq. (3) additionally breaks SU{2)i x SU{2)2 down to the SM SU{2) 
group. 

The hermitian matrix 11 in eq. (6) contains the Goldstone and (pseudo)-Goldstone 
bosons: 

1 



n = 



/ 




V2 


0^ 






H 





H* 






V2 


V2 




\ 





V2 


e 


J 



C°diag(l, 1,-4, 1,1) , 



(8) 



where H — [h^^h^) is the Higgs doublet; is a complex SU{2) triplet given by the 
symmetric 2x2 matrix: 



= 



(9) 



the field C° is a singlet and finally, ^ is the real triplet of Goldstone bosons associated 

to SU{2)i X SU{2)2 SU{2) breaking: 



L 2? 



(10) 



The kinetic part of the Lagrangian is 



A = ^IV[P.S)P'^E)t] (11) 

where ^ 

L>^E = a^E - i gjW^iQ";^ + sgf ) - i^'Sy(yE + Ey^). (12) 

In this model additional fermions are introduced as a vector-like coloured pair 
t^f^ to cancel the quadratic divergence from top loops (we neglect the other Yukawa 
couplings). The relevant part of the Lagrangian containing the top Yukawa coupling 
is given by 

Cf — -XifeijkexyXi^jx^kyu'^ + ^2/^^ + h.c, (13) 

where Xi = {b,t,t), indices i,j,k run from 1 to 3 and x,y from 4 to 5, and e is the 
completely antisymmetric tensor. 

As explained in [2] considering gauge and fermion loops, one sees that the La- 
grangian should also include gauge invariant terms of the form, 

-Cs^V = 2cfg^J2mQtmQtm + 2c/V'Tr[(FE)(FE)*] 

a 

+ 4c7^A?e"'^e,,E,^E,-,E^^*E^'- , (14) 



with c and c' assumed to be constants of 0{1). 

This Lagrangian produces a mass of order / for the gauge bosons (1^') associated to 
the broken (axial) SU{2), for a vector-hke combination of i and u'^ and for the complex 
scalar 0. The singlet C° is a pure Goldstone [associated to the breaking of the U{1) 
symmetry left ungauged] that will play no significant role in the discussion (it can be 
given a small mass to avoid phenomenological problems by adding explicit breaking 
terms). Finally, the Higgs boson gets a small tree level mass of order g'^ and a quartic 
coupling (not suppressed by g'). 
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Figure 1: Periodic structure of the Little Higgs effective potential for T = 0. 

3 Effective Potential (T=0) 

From the Lagrangian (14) we can extract the tree- level effective potential for the real 
scalar field h = ^/2Tie{H). In most previous papers an expansion in h is performed. 
Here we avoid doing this, except to illustrate a few aspects, and keep the full dependence 
on h. The tree-level potential is a function of 



(E)=exp[2z(n)//] Eo 



m\ m\2 
h + iV2^ sin 2a - 4^t4- sin^ a 



h 



(15) 



where (11) is only nonzero through h and a = h/{\/2f). This is a periodic function 
and result the potential 



V{h) = Constant + f 



2cg"sl + 



1 



(16) 



[where A-t = c{gf ± g^) ± 16c'A^] is invariant under h ^ h + mcf. The minimum at the 
origin h = is replicated at invr/ (see figure 1) with barriers of height (2c(7'^ + A+/2)/^ 
separating these minima. In each of them, in spite of appearances (the fact that h ^ 0) 
the electroweak symmetry is unbroken. In fact one can show that the mass spectrum 
is the same in all these vacua, in particular SM gauge bosons are massless. In other 
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words, the order parameter for electroweak symmetry breaking is Mw rather than (h).^ 
Regarding the local properties of these minima, an expansion in powers of h around 
h = gives: 

V{h) = Constant + cg'^fh^ + ^{^K ' ^cg'^)h^ + ... (17) 

As expected, the only^ mass term for h is due to the U{1)y gauge couphng g' . This 
offers an immediate possibility for electroweak symmetry breaking if one chooses c < 0. 
In that case the minimum is at 

{sD = ^ , (18) 

or 

^ ^ - (3A+ - Acg'^) ■ ^''^ 
However, this trcc-lcvcl breaking is problematic. The expansion of the potential to 
fourth order in h and t = \/2Iie{—i(p) is: 

V{h, t) = Constant + cg'^fh^ + ^(3A+ - 4cg'^)h* + -^X.fhH 

+ (A+ + Acg")ft' - ^(4A+ + 17cg")t'h' - ^(A+ + ^cg'Y + - (20) 

In the presence of the coupling A_, a vev for h induces a tadpole for t so that the 
previous minimum gets slightly displaced. Prom (20) we get 

^'^-2^/2(A+ + 4c^'^) P ' ^''^ 

but one cannot obtain | A_ | <^ A+ + 4:cg'^ and as a result {t) turns out to be too large. 

Therefore the breaking has to be triggered by one-loop radiative corrections. In 

order to compute the one-loop potential, both at T = and at finite T, we need the 

mass spectrum in an arbitrary Higgs background h. This is easy to compute and one 

finds that the masses inherit the periodicity of (E) (we give the details in Appendix A) . 

Then the one-loop potential is also a periodic function of h, as shown by figure 1, where 

^The periodicity extends in fact to the plane {Re{h), Im{h)} as the potential is a function of \h\ 
only. The minima at mrf are then circles in that plane, which are nevertheless equivalent to the point 
at the origin h = 0. One can think of Re{h) and Im{h) as orthogonal coordinates on a sphere to 
understand the periodic structure of the potential: different minima correspond to the same point on 
the sphere. 

"^There are additional two-loop contributions to the mass term [5]. 
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Figure 2: Close-up of fig. f to reveal the electroweak minimum. 

it can be compared with the tree level one. In particular, the height of the barriers 
between minima is somewhat reduced. More importantly, the electroweak symmetry 
is now broken: the appearance of symmetry breaking minima is clearer in the close- 
up shown in figure 2. Analytically, this is understood as the result of the negative 
contribution from the heavy top to the Higgs mass. To order h"^ the T = one-loop 
potential reads 



6iV 



647r2 



-12A^M| I log 



A^ 



9 2. .2 A 



g'M' log 



1 



+6 \^X+ + 5cg' 
+6 ( -A , 



A2 



5 ,2 



+ 



(22) 



where 



M| ^ (A? + XI) f , 



/2. r2 



Mt = 2{X+ + 4cg'')f 



Ml 



(23) 
(24) 
(25) 
(26) 
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We have chosen the following values for the parameters of the model: Ai = A2 = 1, 
g-^= g2 = A+ = 0.85 and A_ = -0.0077, so that (h) = 246 GcV for / = 1 TeV.^ 
These parameters also give adequate values for the masses of the non SM particles in 
the model. 

Before proceeding, it is convenient to say a word about our way of computing the 
masses of scalar fluctuations {i.e. the degrees of freedom in E itself). In principle we 
could simply shift hP ^ hP -\- h/y/2 in E to compute /i-dependent scalar masses. To 
do this properly one should take into account that in general, after this shifting, the 
scalar kinetic terms from (11) are not canonical. Therefore one should rescale the fields 
to get the kinetic terms back to canonical form and this rescahng affects the 0{h?) 
contributions to scalar masses. Instead of following this standard procedure we find 
it convenient to use an alternative method which simplifies the calculations and has 
appealing features when one is concerned about the global structure of the potential, 
see below. 

The idea is to treat the new background with {H) — h/\/2 as a basis change (recall 
the discussion of the change from ($) = to ($) = Eq in section 2). The SU{5) 
transformation is now Uh = exp(i(n)//) with Eft = (E) = Uu^qUI- To parametrize 
the scalar fluctuations around this background we again use the exponentials of broken 
generators, but taking into account the effect of the change of basis, which acts on gen- 
erators as Xa UhXaUl- That is, instead of using exp(in//) we use Uh exp{ill/ f )Ul, 
and write for S: 

E = (C/fte'n/^C/J)(C/ftEoC/ft^)(C/*e^n"/^C/ft^) = e^WZ/g^^nz/giW/Zj.^ _ (27) 

Alternatively, defining IIjv = T^a,NX°''^ with Xa,N — exp(i(n)//)Xaexp(— i(n)^//) we 
have 

E = e^"^/^Efee-^"?^/^ = e'^^^/^E^ . (28) 

[This parametrization allows the second equality of eq. (28) to hold as the condition 
Xa,N^h — ^h^a,N ^^^^^ Valid.] The prescription in eq. (27) [or (28)] is to be compared 
with the standard procedure 

E = e2^("+<n»/^Eo . (29) 

^The smallness of A_ gives a small vev to t, which we then neglect to focus only on the /i-direction. 
For the behaviour of the potential along the f-direction see Appendix B. 
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Some comments are in order. First, it is easy to check that with the prescription 
of cq. (27) scalar fluctuations are automatically canonical. Second, we are free to 
choose this parametrization of scalar fluctuations: a general theorem [15] guarantees 
that this different parametrization docs not change the physics. Finally, using eq. (27) 
has another advantage when looking at global properties of the scalar potential. As 
mentioned above, the mass spectrum is the same in all the periodic degenerate minima 
of the potential. For the scalar sector this is trivial to show with our parametrization 
and a bit more cumbersome using the conventional parametrization of eq. (29) which, 
as explained, requires field redefinitions in order to recover canonical kinetic terms. 
Some of these field redefinitions are in fact singular (because some kinetic terms go to 
zero in these minima). Although the singularities cancel out when computing scalar 
masses (because some masses also go to zero in these minima) they can make some of 
the scalar couplings blow up. This can be interpreted as a failure of the parametrization 
(29) to cover the whole parameter space with a single coordinate patch. In this sense 
the parametrization (27) is to be preferred when discussing global properties of the 
scalar potential. 

4 Effective Potential {T^Q) and discussion 

We calculate the finite temperature effective potential at one loop in order to study the 
dynamics of the scalar fields at finite temperature including the interactions with the 
thermal bath. To do this, we include the contributions from all particles that receive 
a correction to their mass from the vacuum expectation value of the nonlinear sigma 
model field (E) = E/^, just as we did for the one-loop potential at T = 0. The one-loop 
thermal integrals for the contributions of bosonic and fermionic particles arc standard, 
sec e.g. [16]. For certain regions of h some scalars might have negative masses squared, 
rn^ih) < 0. For these we simply take the real part of the corresponding thermal 
integral^. 

Figure 3 shows the behaviour of the global structure of the potential for increasing 

values of the temperature. We see that the barrier between two minima decreases 

when T is increased, eventually turning into a local minimum (where the electroweak 

®In [17] an alternative treatment of these integrals, with an infrared momentum cutoff |m^|, is 
used. We have checked that both prescriptions are very close numerically. 
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symmetry is broken) that becomes the global minimum if T is even higher. We can 
define a critical temperature Ti at which this local minimum is degenerate with the 
minima at /i = 0. Numerically we obtain Ti ~ 0.96/ for our particular choice of pa- 
rameters (Ti ~ / will be generic). An analysis of the thermal behaviour of models with 
pseudo-Goldstone bosons was done in [18], in a different context. We find that some 
of the generic features discussed in those papers resemble the behaviour we encounter 
for Little Higgs models. 

In order to understand analytically the particular behaviour of the potential shown 
in figure 3 it is enough to consider the case olT > f. At such high T (still below A) 
an expansion in powers oim/T gives a good approximation to the potential. Writing 
also the one-loop T = part, one obtains in general 

^2 



24 V " ' 2 " "^J 12% "' ^ 



1 

+ 



647r2 



IVM^ [log ^ + kb)- IVM) flog ^ + KF 



(30) 



where B — {S, V} labels bosonic degrees of freedom (scalar and vector bosons) and 
F fermionic ones. The constants k,b,f come from the one- loop T — potential and 
are simply ks — kf — —3/2 and Ky — —5/6 while cb,f come instead from the finite 
temperature part and are logce — 5.41 and logcp — 2.64. Keeping only the dominant 
term we get for our model (see appendix A): 

StV = -^^fT' { [l2A+ + A?(64c' + 6) - ^igl + g^)] 4 + g" [48c - | «^ } ■ (31) 

First we see that expanding around h = 0, the thermal mass is proportional to g', 
as expected on general grounds. For the range of temperatures we are considering, 
this thermal effect is not strong enough to produce a minimum in h — 0, the one- loop 
T — corrections dominate and indeed the potential at h — has a maximum. Next 
we examine larger values of h. In eq. (31) the gauge contribution is subdominant, 
and the scalar and fermionic terms provide a deep minimum for the potential at the 
intermediate point a — 7r/2. This is precisely the high T behaviour that we find in 
figure 3. 

The results indicating symmetry nonrestoration of the gauge electroweak symmetry 
directly lead us to pose the question of what could be the associated cosmology. As we 
have said already, we are restricted by our analysis to temperatures beneath A/tt ~ 4/. 



12 



Suppose now that we follow the behaviour of the system starting from low temperatures 
and heating up the thermal bath. As one can infer from the lower plot of Fig. 3, the 
electroweak gauge symmetry is restored at some temperature Tc„ of the order of the 
electroweak scale. This means that h = (modulo tt/) becomes the ground state of the 
system. This is completely analogous to what happens in the SM. However, a further 
increase of the temperature results in a decrease of the height of the barrier between 
the periodic minima at h — 0, tt/, .... For even larger values of the temperature beyond 
7i ~ /, the maximum at h — irf /2 (modulo tt/) turns into the global minimum and 
the barrier among the new equivalent vacua increases with temperature. 

The cosmology of the Little Higgs model will therefore depend strongly on the 
maximum temperature that the Universe has attained after a period of inflation [19], 
which is necessary to explain the homogeneity and the isotropy of our observed Uni- 
verse. Suppose that such a temperature is smaller than Ti and that the minima at 
h = 0,71 f, ... are always the ground state until the electroweak phase transition. If so, 
the situation is completely analogous to the SM. Suppose, however, that the highest 
temperature after inflation is indeed larger than Ti. Under these circumstances, the 
Universe is characterized by a high temperature phase during which the electroweak 
symmetry is broken. This period will last till the temperature drops somewhat below 
Ti at which point there will be a phase transition to the symmetric phase (followed 
later on by the usual SM electroweak phase transition back to the broken phase) . This 
prediction is completely different from what is obtained in the SM. It would be very in- 
teresting to analyze the role played by this non-standard phase of broken symmetry in 
the evolution of the early Universe, e.g. on the possible baryon asymmetry production 
[20]. 

In view of the small thermal mass that comes from (31) (which is even zero in most 
Little Higgs models), one might wonder if two-loop T 7^ effects can then become 
important^. In fact, at T = 0, two-loop quadratically divergent corrections to the 
Higgs mass parameter go as Sm"^ ~ A^/(167r^)^ ~ /^/(IGtt^), the same order as one- 
loop non-divergent contributions. (We did not include them because the one-loop 
contributions will in principle dominate due to logarithmic enhancement). At T 7^ 
one would correspondingly expect corrections of order ^fn^ff ~ ^^/(ISvr^), which may 
be important precisely when the one loop thermal mass is suppressed (or absent). In 
^We thank Bob McElrath for pointing this out. 



13 



fact, as the T = Higgs mass squared is itself of order p j (167r^), a two-loop thermal 
mass of that order can be relevant when T > f . A complete two-loop calculation of the 
finite T effective potential is beyond the scope of this paper, but it is of interest and we 
plan to undertake it in a future analysis. At this point we simply make two remarks. 
The first is that, for values of h away from the origin, the one- loop T"^ contribution 
to the potential [eq. (31)] is negative and not suppressed. Therefore, the potential 
in that region will not change much after including two-loop corrections, which will 
be sub-dominant there. This gives us confidence on the inverse symmetry breaking 
behaviour we have found. The second remark is that we see no reason to expect that 
the two-loop contributions to the Higgs thermal mass will be positive. If they are, the 
details of the transition around Ti will change (but not the existence of the transition 
itself) while if they are negative, our conclusions would be even stronger. 

Let us close this section by a couple of comments. First, a complete study of the 
scalar potential should also include the temperature evolution of the triplet vcv. For 
that purpose the spectrum in a more general background with nonzero h and t is 
needed. This complicates significantly the analysis, especially if one insists on keeping 
the h and t dependence to all orders (the potential is also periodic along the t direc- 
tion). We nevertheless performed such analysis and some of our results are presented 
in appendix B. Schematically, the potential has an egg-crate structure with different 
barrier heights for each direction. Although the potential at finite temperature along 
the t direction can behave quite differently, depending on the choice of parameters, 
from what we described for the h direction we find it interesting at present to focus on 
the Higgs direction alone. 

Secondly, the previous analysis has assumed that / is a constant, independent of 
temperature. However, we know that / is in fact the vev of E (along a particular 
direction), producing the breaking SU{5) — > SO{5). As such it is a dynamical variable 
and it is expected that at sufficiently high T one will get / — > 0, corresponding to a 
critical temperature for the SU{5) — > S0{5) transition. As we do not know what is 
the physics beyond the cutoff scale A, where the dynamics of the system is superseded 
by the UV completion of the theory, we do not know the potential that produced the 
vev / in the first place. Therefore, it is difficult to be precise about T*. However we 
can make an estimate of the temperature behaviour of / (in the spirit of [17] for the 
more complicated case of the chiral condensate in QCD). We can just approximate the 
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zero temperature potential for / by a Mexican hat potential 

v{f) ^ K'if - {frr , (32) 

where k is some unknown constant and (/) is assumed to be ~ 1 TeV (it corresponds 
to the / used in the rest of the paper). We then add to the potential (32) the finite 
temperature corrections coming from all the particles that have an /-dependent mass. 
These are hsted in eqs. (24)-(26). The value of k, in eq. (32) is now crucial to the change 
of / with increasing T. We argue that a natural choice is k ~ 27r, because fluctuations 
around the minimum of the potential (32) along the / direction have mass 2fi;(/) and 
we are assuming that the only scalar fields below the cutoff scale A = 47r(/) are those 
contained in 11. Therefore we should demand 2«;(/) > A which indeed translates into 
K > 27r. Choosing then that value of k for the numerical evaluation of the effective 
potential for /, we find that / changes very little with T. For the extreme value 
T = 4(/) we obtain /(T) ~ 0.9/, just a 10% decrease. 

5 Conclusions 

We have shown, on very general grounds, that the behaviour of Little Higgs models 
at finite temperature is considerably richer than in the Standard Model. In particular 
we have studied the effective potential at finite temperature which, as the Higgs h in 
these models is a pseudo-Goldstone boson, is a periodic function of h. 

Although the electroweak phase transition is expected to occur just like in the 
Standard Model, at higher temperatures, when the new states introduced in Little 
Higgs models at a scale / ~ 1 TcV arc thermally produced in the plasma, the his- 
tory of the early Universe changes dramatically. At some temperature Ti ~ / a new 
minimum where the electroweak symmetry is broken becomes the global minimum of 
the potential and the gauge electroweak symmetry becomes more and more broken as 
the temperature continues to increase. Being a perturbative statement, it would be 
interesting to see if such a behaviour persists when non-perturbative corrections are 
accounted for, or what is the fate of this broken minimum in the context of the UV 
completion of the theory. At any rate, this rich structure at high temperatures might 
have cosmological implications which arc worth studying. 

As a spin-off, we have introduced a new parametrization of scalar fluctuations 
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around displaced vacua in nonlinear sigma models which has very appcaUng features 
and might be of interest for phenomenological studies. 

Acknowledgments 

We acknowledge very useful discussions with (and help from) Alberto Casas, Mikko 
Laine, Bob McElrath, Jesus Moreno and Ann Nelson. This work was supported in 
part by Colciencias under contract no. 1233-05-13691. J.R.E. thanks CERN and 
Univ. of Padova for hospitality during several stages of this work. M.L. thanks the 
I.E.M. of CSIC-Madrid and LPT-Universite de Paris Xl-Orsay for hospitality during 
the completion of this work. A.R. thanks CERN for hospitality. 

A. Spectrum of the Little Higgs model 

In this appendix we present the mass matrices in a Higgs background for the Little 
Higgs model [10] studied in this paper. These masses are needed in the calculation of 
the one-loop Higgs potential, both at zero and finite temperature. We keep the exact 
dependence on the Higgs background = h/\/2 so that we can study the global 
structure of the potential. We also present the mass eigenvalues in an expansion up 
to . The scalar sector contains a complex SU{2) triplet, a real triplet, the Higgs 
doublet and a singlet, a total of 14 degrees of freedom that are distributed in one 
doubly charged scalar, three charged fields and 6 neutral and real fields (4 CP odd 
and 2 CP even) . The gauge boson sector contains a heavy W and a heavy Z' besides 
the SM gauge bosons. The fermion sector of relevance for our calculation contains 
a heavy top in addition to the SM top quark. We follow the notation introduced in 
section 2. We remind the reader here that A+ = cG\ + 16c' Af and A_ = cG'i, — 16c' Af, 
with = gj + = gj — 92- We also use a = h/{\/2f), where h is normalized 

as a real field. (At T = one has {h) = 246 GeV.) The mass matrices presented 
below correspond to canonically normalized fields (see section 3 for a discussion on this 
point for the scalar sector). The mass eigenvalues of all scalar fields are periodic under 
q; — >■ a -|- TT. 
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Doubly charged scalar 

The field 0"*"+ has a mass 

Mj+, = 2/2(A+ + 4c^'^ - IQdXlsi) = 2(A+ + Acg")f + 0{h'/f) . (A.l) 
At order h? this state does not contribute to the trace of the mass squared operator. 

Charged scalars 

In the basis {^"'", 0+, /i"*"}, the mass matrix for the charged scalars is 



-2A_|_Cq,(2 + Ca)s'^i2 



2A+Ca(2 — Ca)c^i2 ~;^A_S2aC; 



,2 

•a/2 



^A_S2a4/2 



A q2 2 



For the trace we obtain 



2(1 - 3c„)s^/, 4 
si 2(l + 3cJc2/2 











2ci 



(A.2) 



2IVM^ = 4A+/^(1 - 4) + 4c^'V'(5 - 4.^) , 
and an expansion in powers of h gives 

2TrM| = Constant - Scg''^h^ + 0{h'^ / f) , 



(A.3) 



(A.4) 



with the only contribution to order h'^ being that of the U{1)y sector. A similar 
expansion for the mass eigenvalues gives 



ft 2 
mJ, = 2(A+ + 4c/)/^-i 



A+ + 5cg'^ - ^ 



,/2 



= 2cg''f + ^ 



A+ - 2cg'' - 



A+ + ?>cg' 



A+ + ?>cg' 



h' + ... 



,/2 



(A.5) 
(A.6) 

(A.7) 



For h = 0, = is the charged Goldstone boson associated to the gauge symmetry 
breaking SU{2)i x SU{2)2 SU{2)sm and we get = 0. 
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Neutral scalars 



For the neutral scalar fields </7° we use the basis {('^ , ^'^ , (f)'^^ , (f/''^ , h'^^ , hP'^} , where we 
have decomposed the complex fields 0° and hP as /i° = {h^^ + ihP^) / V2 and 0° = 
+i(f)^) I -\/2. The mass matrix for neutral scalar fields, Mq , breaks up in two blocks: 
one for the pseudoscalars ^ ■, 0°*, and the other for the scalars {0°'', The 
4x4 block is 



^1x4 



eg f 



10X^2 -sl)sl 
V2X^{2-sl)sl 








10A+(2 - sl)sl V2A+(2 - sl)sl -2X+{Acl - 4) 4A_c„S2a 





54 V^sl 



4A_CaS2o 



— A+s 



lOs 







TOsl V2sl -2{A - 5sl) 
-2d 



V2sl 



(A.8) 



which has a zero eigenvalue, corresponding to the eigenvector (1, — v^S, 0, 0)/\/6 [the 
neutral Goldstone boson associated to the breaking of an ungauged ^/(l)]. 
The 2 X 2 block is: 



^2X2 



f 



X42cl-si) X43sl-2)s^ 
A_(34-2)s« A+(3-44)4 



+ 2cg'y 



4 - 



C2a 







(A.9) 



Prom these matrices we obtain 



TrM-' = 4A+/^(1 - 24) + Acg"f{5 - 8^) . 



(A.IO) 



An expansion in powers of h gives 

TrM^ = Constant - 16cg'^h'^ + Oih'^/f) , 



(A.11) 



with the only contribution of order h'^ being that of the U{1)y sector. The expansion 
of the mass eigenvalues gives, for the 4x4 block: 



<o = 



m^o=2(A+ + 4c/)f - 



A+ + hcg'"^ - 



A2 



A+ + 2>cg' 



/2 



2cg''f + ^ 



A 



+ 



2cg 



l2 



2A2 



A+ + Scg' 



/2 



(A.12) 
(A.13) 

(A. 14) 
(A.15) 
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while for the 2x2 block we get 

mJo = 2(A+ + 4c/)/^- 



A+ + hcg'"^ - 



mlo = 2cg"f + - 



3A+ - 4c^ 



/2 



A+ + 2>cg'^ 
2A2 



X+ + 3cg'^_ 

For h — we have two zero mass eigenvalues: the Goldstones and 



(A.16) 
(A.17) 



Gauge Bosons 

The mass matrix for the charged gauge bosons {W^, W'^) in the interaction basis is 
given by 



gl -9i92cl 

r3. 



gl 



-gig2c^ 

which has periodic mass eigenvalues. 

There are three neutral gauge bosons: Z^, Z'^ and the photon. Their 3x3 mass 
matrix in the interaction basis is 



f2 

ZA g 



gl{2 + si) -gMl + ct) 2g,g'sl 
-g,g,{l + ci) gi{2 + si) 2g2g'sl 



2g,g'sl 2g,g'sl Ag'^sl 

again with periodic eigenvalues. It can be easily checked that annihilates the 

vector {1/ gii l/5'2, —^Ig'), which corresponds to the photon. For the trace of we 
get 

(A.20) 



(A.19) 
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and expanding in powers of h, 



TrM^ = Constant + -g'^h^ + 0{h^ / f) 



The expansion of the mass eigenvalues is 



m^± = 








m|o = 









(A.21) 

(A.22) 
(A.23) 
(A.24) 
(A.25) 



where we have used g = glgl/{gl + gl) 
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Fermions 

For the top quark and its heavy partner we have 

which has periodic eigenvalues. The fermionic trace is 

T,{MfM}) = f[\l + \l{l-sl)\ , (A.27) 

which does not contribute to order h?. The expansion of the two mass eigenvalues can 
be written as 

ml = h^h' + ... (A.28) 
ml = {XI + XI) f - \xlh\.. (A.29) 

where the top Yukawa is given hy X^ = 2X\X\/{X\ + X\). 

B. The triplet direction 

As mentioned in the main text, in this model a nonzero vacuum expectation value of 
{h^) = h/ \/2 will induce a vev for the triplet field (0q) = t/ The tree-level potential 
to all orders in h and t is given by, 



4 



2 



(t2 + h^) 

+ h'f'sls,, + \thl, + h^4] , (B.l) 



(t2 + /,2)2 V y ■ (^2 + /,2) 

where P — {t'^+h^y^^/ (V^f). We see that the potential is periodic along the ^-direction 
as well. 

A complete analysis of the one- loop corrections to the potential (B.l), both at T = 
and T 7^ requires the calculation of the mass matrices of all particles in the presence 
of the background fields h, t, which is not an easy task, although we did perform it. The 
generic expressions one obtains for these mass matrices are too long and complicated 
to be given here. 

In this paper we limit ourselves to give an example of the kind of dynamics one 
might obtain along the triplet direction. The interest of this is not hmited by the fact 
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that the phenomenological comphcations associated to the presence of a triplet vacuum 
expectation value can be neatly solved by imposing a T-parity on the model [11]. As 
we will explain next, even if one does not have a tadpole for t the triplet direction 
might play a role in the thermal evolution of the Universe. 

The argument that leads one to expect a peculiar behaviour of the effective po- 
tential at finite temperature along the Higgs doublet direction (namely the fact that 
quadratically divergent corrections to the Higgs mass vanish) does not hold for the 
triplet direction. Explicitly, an expansion of the t-dependent masses for the gauge 
bosons along the triplet direction gives 

m^± = ^gH'' + ... (B.2) 

mU = \Glf - \gH\.. (B.3) 

m|o = (/ + /)t2 + ... (B.4) 

ml^-\Glf + \^/... (B.5) 

so that the trace is 

3TrM^ = Constant + ^{G\ + Ag''^)t^ + ... (B.6) 

For fermions one gets 

12TrM| = Constant - 2AXlt^ + (B.7) 

while for scalars 

TrMi = Constant - -(318ct/'^ + 1488c'A? + 4Q\+)t^ + ... (B.8) 
5 

These results determine that the thermal potential along the triplet direction has 
a contribution 

SV{T) = —nlt^T'^ 



1 

48 



^{Gl + Ag'^) - ^(318c/ + 1488c' + 49 A+) - 2AXl 



t'T' . (B.9) 



The negative contributions to the traces in (B.7) and (B.8) dominate in the potential 
(B.9), making > 0. This unusual behaviour implies that for temperatures above a 
critical value 

= Mt/Kt , (B.IO) 
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where = 2{X+ + 4cg' )p is the triplet mass squared, the electroweak symmetry 
gets broken along the triplet direction. That is, even if a T-parity forbids a tadpole 
h?t in the potential, and even if the mass of the triplet is large (of order /), thermal 
corrections will trigger electroweak symmetry breaking along the triplet direction, for 
some f. Of course one can increase the value of by appropriately choosing the 
parameters of the model and in this spirit we have not discussed this transition in the 
main text, and have focused instead in the Higgs direction. Nevertheless, this makes 
the possible structure of phase transitions in these models even richer. 
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